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Abstract. We discuss mass transfer in eccentric binaries containing a white 
dwarf and a neutron star (WD-NS binaries). We show that such binaries are 
produced from field binaries following a series of mass transfer episodes that 
allow the white dwarf to form before the neutron star. We predict the orbital 
properties of binaries similar to the observed WD-NS binary J1141+6545, and 
show that they will undergo episodic mass transfer from the white dwarf to the 
neutron star. Furthermore, we describe oil-on- water, a two-phase SPH formalism 
that we have developed in order to model mass transfer in such binaries. 



1 Introduction 

Interacting compact binaries are responsible for a wide range of astronomically 
interacting events. Accreting white dwarfs (WDs) are believed to be the pro- 
genitors of type lA supernovae. Cataclysmic variables also occur when white 
dwarfs accrete matter but under different conditions. Higher-mass accretors - 
neutron stars (NSs) or black holes (BHs) form X-ray binaries. A recent paper 
( Metzger 20 111 ) suggests that the tidal disruption of a WD by either a NS or a 



BH may lead to an accretion flow that is dominated by nuclear burning, with a 
potentially unique nucleosynthetic signature. 

Here we consider the problem of modeling the onset of mass transfer in 
an eccentric binary containing a white dwarf and a neutron star, a potential 
progenitor of such a nuclear-dominated accretion flow. We discuss the formation 
of such binaries and show that when they come into contact they are expected 
to undergo episodic mass transfer. Finally we outline the oil-on-water SPH 
technique that we have developed to study these binaries. 



2 Formation of eccentric WD NS binaries 



There are two clear detections of eccentric binaries containin g a white dwarf 
and a neutron star ; B2303-I-46 ( van Kerkwijk &: Kulkarni 19991 ) and J1141-6545 



20001 ). One might expect that such binaries would not exist, as 
the neutron star, having a more massive and hence rapidly evolving progeni- 
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tor, should form first. Tidal circularisation during the giant phase of the white 
dwarf's progenitor would then remove any eccentricity imparted by the super- 
nova. Therefore, to produce such a binary, the white dwarf must form firs t. 



A formation mechanism was outlined bvlPortegies Zwart Sz YungelsonI (11999 ) 
and developed bv lTauris Sz Senne li (|2000l ). lDavies et all (|20oi )" andlChurch et al 



(|2006i ). In it, the binary initially contains two main-sequence stars of similar 
but not identical, mass, and has an orbital separation of around 100 i?0. The 
initially most massive star fills its Roche Lobe in the Hertzsprung gap and trans- 
fers mass stably to its companion. Subsequent phases of mass transfer can follow 
several pathways, but in all cases the initially most massive star fills its Roche 
lobe again, transfers the majority of its mass to its companion and finally be- 
comes a white dwarf. The companion star is now massive enough to evolve 
a Chandrasekhar-mass degenerate core, which collapses leading to a supernova 
and, if the resulting binary is bound, the formation of an eccentric binary. We 
refe r to such systems as WD-NS binaries. 

Church et al.l ( 20061 ) predict the binaries similar to J1141-6545 form when 



a fourth episode of mass transfer occurs after the white dwarf has formed. This 
mass transfer is inevitably from the most massive star at this point (the pro- 
genitor of the neutron star) and hence causes the orbit to shrink, producing a 
tighter binary after the supernova. Such close binaries merge in less than the 
Hubble time and hence are of interest to us. 

Following the formation of the binar y, the emiss ion of gravitational waves 
causes its orbit to shrink and circularise ( Peters! 19641 ) . In the fully circularised 
case, mass transfer from the white dwarf begins once it fills its Roche lobe. How- 
ever, these binaries still have a residual eccentricity at this point. In Figure[T]we 
plot the orbital properties for a sample of WD-NS binaries that for m via a path- 



way similar to that inferred for J1141-6545, as computed by BSE (iHurlev et al 



[2002 ). It can be seen that our simulations produce binaries with orbits similar 
to that of Jl 141-6545. The residual eccentricity when these binaries come into 
contact is typically between 10~^ and 10~^. As the surface pressure scale height 
of white dwarfs is very small these eccentricities are large enough for episodic 
mass transfer to occur. 



3 Using smoothed-particle hydrodynamics in eccentric binaries 



In order to model episodic mass transfer in eccentric WD-NS binaries we have 
developed an extension to SPH which we call oil-on- water. SPH is a Lagrangian 
particle-based method for s olving the equations of hydrodynamics ( Lucvl Il977l: 
iGingold Sz Monaghan 19771 ). The implementation that we use follows Bend 
(| 19901 ). In vanilla SPH, one requires a very high mass resolution to model mass 
transfer between stars. For ex ample, in a cata clysmic variable the mass transfer 
rate is 10"^ to lO-^Moyr"^ (|Pattersonlll984l ): given a typical orbital period of 
1 day the mass transferred per orbit is one part in 10^^ of the total mass. For 
studies of continuous mass transfer in circular binaries this is not a problem, as 
the fiow is steady and hence can be considered over a large number of orbits. In 
the case of episodic mass transfer that switches on and off each orbit, however, 
a new approach is required. 



Mass transfer in WD-NS binaries 
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Figure 1. Left panel: The eccentricfties of our simulated sample of WD- 
NS binaries as a function of their semi-major axis, both at formation time, 
i.e. just after the supernova. The filled triangle represents J1141+6545. Right 
panel: Eccentricities of WD-NS binaries at the point when the white dwarf 
fills its Roche lobe, as a function of the time taken for them to spiral together. 



To avoid this problem we split the star up into two separate phases. Heavy 
"water" particles form the body of the star and account for essentially all of 
its mass. Light "oil" particles form an atmosphere that is transferred during 
interactions between stars. As the mass of the oil particles is much less than 
that of the water particles numerical stability requires them not to come into 
contact. To prevent this from happening we introduce an additional force to 
stop the oil particles sinking into the star. This additional force, acting on an 
oil particle of mass rrii , depends on the water particle number gradient n as 



Fn = fi^^f^Vn; Vn = Y: VW°^{nj, h^j) (1) 



where /? is an adjustable parameter. We utilise a different kernel with a 
break point at a smaller multiple of the smoothing length h in order to keep the 
oil layer closer t o the surface of the star. Full details of the method are given 
in IChurch et al.l ( 20091 ). It should be noted that the oil particles are not test 



particles; they obey the full SPH equations when interacting with one another, 
which lets us follow the processes inside the accretor's Roche lobe. 

To test the validity of the model we investigated mass transfer from a low- 
mass star to a white dwarf in eccentric binaries. For the body of the star we used 
15 390 water particles with a total mass of 0.6 Mq. The atmosphere was made 
up of 39 691 oil particles with a mass of 10"^'* Mq each. This stellar model was 
placed into a binary with a 1 Mq white dwarf, modelled as a point mass, and 
relaxed at apocentre. The orbit was then made eccentric. For each eccentricity 
we varied the semi-major axis in order to find the pericentre separation at which 
mass transfer started. We found that mass was transferred on each orbit shortly 
after pericentre passage, and that the semi-major axis required for mass transfer 
decreased with increasing eccentricity, compared with the value that would be 
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Figure 2. Left panel: structure of an oil-on-water star. Particles within 
h of the X — y plane only are plotted. The oil particles are plotted with 
larger dots. A clear atmosphere can be seen within the oil layer. Right 
panel: the ratio between the periastron separation required for mass transfer 
assuming that the star must fill its Roche lobe at periastron and that measured 
from our simulations, plotted as a function of eccentricity. The minimum 
separation required for mass transfer to take place is seen to decrease linearly 
with increasing eccentricity. 



implied by requiring the star to fill its Ro che lobe at peri c entre (see Figure [5]) . 
This is in agreement with the results of Sepinsky et al. ( 20071 ). who take an 



analytical approach to determining the point of onset of mass transfer. 
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